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[ — ' The interaction of dust grains originating from the local interstellar cloud with the environment 
inside the heliosphere is investigated. As a consequence of this interaction the spatial 
distribution of interstellar dust grains changes with time. Since dust grains are charged in the 

>■ ■ interplanetary plasma and radiation environment, the interaction of small grains with the 
S ■ heliosphere is dominated by their coupling to the solar wind magnetic field. The change of the 
ff^ , field polarity with the solar cycle imposes a temporal variation of the spatial distribution and the 
^ I flux of small (radius smaller than 0.4 /im) interstellar dust grains in the Solar System, whereas 
Qv ' the flux of large grains is constant because of their negligible coupling to the solar wind magnetic 
Q\ ■ field. The flux variation observed by in-situ measurements of the Galileo and Ulysses spacecraft 

r-| ! are reproduced by simulating the interaction of interstellar grains with charge-to-mass ratios 

Q-t\ between 0.5 C kg~^ and 1.4 C kg~^ with the interplanetary environment. 
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1. Introduction 



2. Model Description 



Levy and Jokipii [1976] showed that the dominant 
force on small dust grains in the Solar System is the 
Lorentz-force caused by the solar wind magnetic field 
sweeping by the dust grains, which are charged to a 
positive surface potential, mainly by electron emis- 
sion due to solar UV photo-effect. It was concluded 
that due to the coupling to the radially outward flow- 
ing solar wind, small interstellar grains do not exist 
inside the heliosphere which defines the domain of the 
solar wind. For big grains, Lorentz-force is less im- 
portant, because they have smaller charge-to-mass ra- 
tios, but from the determination of the size range and 
size distribution of interstellar grains in our galaxy 
by fitting the wavelength dependence of extinction of 
starlight on its way through the interstellar medium 
[e.g. Mathis et al, 1977], it was beheved that typi- 
cal interstellar grains are not larger than 0.5 ^m in 
diameter. So the conclusion was drawn that there 
are no big interstellar dust grains and the small ones 
are swept out of the Solar System by the solar wind 
magnetic field. Gustafson and Misconi [1979] on the 
other hand, argued that the electromagnetic inter- 
action enhances the spatial density of small (radii 
between 0.0005 iiva and 0.12 fim) interstellar grains 
upstream of the Sun. The detection of interstellar 
grains in the Solar System by the dust detector on- 
board Ulysses [Griin et al., 1993] and its confirma- 
tion by data from the Galileo dust detector [Baguhl 
et al., 1995] proved the existence of grains larger than 
0.5 fim and that smaller grains can penetrate into 
the Solar System, although they are depleted in num- 
ber [Griin et at, 1994]. It was concluded [Griln et 
al., 1994; Gustafson and Lederer, 1996; Grogan et al., 
1996] that the Lorentz-force on interstellar grains not 
only repels but also focuses them, depending on the 
configuration of the solar wind magnetic field, and 
thus on the phase of the solar cycle. The objective 
of this work is to quantitatively determine the spa- 
tial distribution of interstellar dust grains in the So- 
lar System as a function of the solar cycle. To meet 
this objective, the motion of interstellar grains in the 
Solar System is simulated by numerically integrating 
the equation of motion for a big {N — 10^ to 10^) set 
of interstellar grains. In section 2 the assumptions 
and approximations used in the model are described. 
Section 3 presents the resulting three-dimensional and 
time-dependent spatial distributions of grains in the 
Solar System. The results are discussed and com- 
pared with the in-situ measurements of the Ulysses 
and Galileo spacecraft in section 4. 



As initial conditions, a homogeneous, mono-direc- 
tional stream of grains is assumed arriving from an 
upstream direction of Aoci — 259° (heliocentric eclip- 
tic longitude) and j3cc\ — +8° (heliocentric ecliptic 
latitude). This direction is equal to the upstream di- 
rection of interstellar dust grains in the Solar Sys- 
tem as it was determined by analyzing the direc- 
tional information of the Ulysses in-situ measure- 
ments [Baguhl et al, 1995; Frisch et al., 1999]. This 
direction is also compatible with the upstream direc- 
tion of interstellar neutral helium atoms, determined 
by the measurements of the Ulysses/GAS experiment 
[Witte et al, 1996]. It was shown [Griin et al, 1994], 
that the impact velocity of interstellar grains mea- 
sured by Ulysses after Jupiter fly-by was compatible 
with the velocity of the interstellar helium atoms. Be- 
cause the impact velocity of an individual grain is de- 
termined with an uncertainty of factor 2 by the in-situ 
measurement [Griin et al, 1992], the helium velocity 
of Woo = 25.3 km s~^ was adopted for the model. In 
the simulation, a "wall" of grains is placed at a he- 
liocentric distance of 100 AU every time interval Ai. 
The wall covers an area of 30 AU x 30 AU for big 
and 120 AU x 200 AU for small grains, and is v^o^t 
thick. Initially the grains are distributed uniformly 
inside the wall. 

The position x and velocity x of each grain as 
a function of time is determined by integrating the 
equation of motion given by: 
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where (3 is the strength of radiation pressure expressed 
as the ratio of the magnitude of radiation pressure 
force to the magnitude of gravity, G is the gravita- 
tional coupling constant, Mq the mass of the Sun, 
q and m the charge and mass of the grain, Wsw the 
solar wind velocity, and _Bp the solar wind magnetic 
field. In the case q/m -^ 0, i.e. if the Lorentz-force 
on the grains is negligible, the spatial distribution of 
grains can be determined analytically [Landgraf and 
Miiller, 1998] and the solution is rotationally sym- 
metric about the axis parallel to the upstream di- 
rection. If the Lorentz-force can not be neglected, 
the solution to equation (1) has no continuous sym- 
metry and the dust grain trajectories are essentially 
three-dimensional. Equation (1) contains three ba- 
sic parameters: the strength of radiation pressure /3, 
the charge-to-mass ratio q/m, and the strength of 
the radial i?r,o and azimuthal B^,o components of 



the solar wind magnetic field. To reduce the num- 
ber of free parameters in the model, spherical grains 
with a constant density of p = 2.5 g cm"^ are as- 
sumed. In this case, /3 and q/m can be expressed 
as a function of grain radius a. In the model the 
calculation of /3 as a function of grain radius by 
Gustafson [1994] is used (see also table 5 below). In 
this Mie-theory calculation, the optical constants of 
"astronomical silicates" introduced by Draine and Lee 
[1984] are assumed. The parameter /3 reaches its max- 
imum for a = 0.18 /im and is greater than unity for 
0.09 /im < a < 0.4 /im. In this regime, interstellar 
grains are effectively repelled by the Sun. For a ^ 0, 
(3 reaches a constant value of w 0.1, and for a — > oo, 
13 depends on a like a^^, which is valid in the limit of 
geometric scattering. 

For the charge on a dust grain in interplanetary 
space it is assumed that the charging process due to 
photo-emission of electrons induced by solar UV pho- 
tons partially balanced by the interplanetary plasma 
environment produces a positive surface potential of 
U = +5 V on the grain surface [Griin et ai, 1994]. 
The charge of the grain is then given by 

q = AireoUa, 

where eo is the permeability of the vacuum. Thus, 
q/m oc a~^, which means that electromagnetic in- 
teraction with the solar wind magnetic field is only 
important for small grains. 

In the model, Parker's [1958] model is used, which 
has been found to be a good large scale approxi- 
mation of the observed magnetic field in interplan- 
etary space [Mariani and Neuhauer, 1990; Burlaga 
and Ness, 1993; Smith et al., 1995]. The radial, az- 
imuthal, and normal components (i?^, S^, Bq) of the 
Parker-field can be expressed as a function of location 
in heliographic coordinates (rhg, Ahg, /3hg): 
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3 nT are the radial and az- 



imuthal field strength at rg = 1 AU, respectively 
[markcitcGMsto/son, 1994]. The Lorentz-force Fi^ on 
a dust grain in interplanetary space caused by the 
solar wind magnetic field is dominated by the ef- 
fect of the field sweeping by the grain with the so- 
lar wind speed, because the grain velocity is much 



smaller than the solar wind velocity. Therefore, the 
Lorentz-force is governed by the term (q/m){B^Vsw) 
(see equation (1)). Although it would seem that the 
Lorentz-force depends on the solar wind speed, this is 
not the case, because at higher solar wind speeds, the 
Parker spiral becomes steeper, and consequently the 
azimuthal component B^ is reduced [Gustafson and 
Misconi, 1979], effectively keeping the Lorentz-force 
constant. The observed sector structure of the po- 
larity of the solar wind magnetic field has been mod- 
eled by Alfven [1977] using a "ballerina model" . This 
model describes a warped current sheet that separates 
opposite polarities of the magnetic field in interplan- 
etary space. In the ecliptic, opposite polarities are 
observed in separate sectors, the boundaries of which 
are given by the intersections of the current sheet with 
the ecliptic plane. For dust grains in the considered 
size-range between a — 0.1 /im and a = 0.7 /im, the 
instantaneous polarity can be replaced by an aver- 
age polarity p which is a function of the heliographic 
latitude /3hg- This averaging process is a good approx- 
imation for grains for which the Larmor-frcqucncy is 
orders of magnitude smaller than the solar rotation 
frequency loq, therefore. 
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where Bq = 3 nT is the interplanetary magnetic field 
strength at 1 AU and ujq = 2.7 • 10^® is the solar 
rotation frequency (rotational period of 27 days). The 
effect of the sector structure on very small (6 nm to 
10 nm) dust grains which have been ejected by the 
Jovian system, has been observed by Ulysses [Zook et 
at, 1996] . The average polarity p seen by a dust grain 
at heliographic latitude /3hg, is given by 



2 . tan/3i 

— arcsm 

TT tane 



hg 



(3) 



Br 



Br,<p — P\B,, 



r,if>\ 



where e is the tilt angle of the (equivalent) flat cur- 
rent sheet, that rotates with the 22-year period of the 
solar cycle [Gustafson, 1994]. From equation (3) it is 
clear that the average field strength Br,^, is small at 
low heliographic latitudes, whereas the instantaneous 
field strength is high at low heliographic latitudes (see 
equation (2)). This is because a particle at low lati- 
tudes experiences a negative polarity nearly as long as 
a positive polarity. It was shown [Morfill and Griin, 
1979] that the stochastic perturbation of the trajec- 
tories of small dust grains by the sectored solar wind 



magnetic field produces an outward directed diffusion 
of these grains. This is not taken into account by the 
model, because the sector structure was removed by 
using an average polarity. 

The average polarity p defines the sign of the az- 
imuthal component of the average magnetic field and 
thus the sign of the Lorentz-force, and because the 
average polarities in the northern and southern hemi- 
spheres are always opposite, the Lorentz-force acceler- 
ates dust grains away from the solar equatorial plane 
for p > and towards this plane for p < 0. The se- 
quence of this focusing/defocusing cycle is shown in 
table 1. 

To simulate the motion of interstellar grains in the 
heliosphere, equation (1) was integrated by using a 
simple Runge-Kutta integrator of fourth order [A^m- 
merical Recipes, 1992]. After the first grains have 
completely traversed the heliosphere, the spatial dis- 
tribution of dust grains has been determined at each 
time step by counting grains inside each cell of an 
orthogonal, equidistant grid. The resulting spatial 
distributions of interstellar grains in the heliosphere 
are discussed in the following sections. 

The model described above uses some non-trivial 
approximations. For the determination of the basic 
parameters /3 and q/m, a spherical grain shape has 
been assumed. On the other hand, the observation of 
the polarization of starlight (for summary see Li and 
Greenberg [1997]) has shown that interstellar dust 
grains (at least the part of the population that causes 
the extinction) must have elongated shapes. A devi- 
ation from a spherical shape would increase q/m for 
the same surface potential U, because the charge car- 
riers on the grains are separated by larger distances. 
The /3- value would also increase for non-spherical par- 
ticles, because it is proportional to the cross-section- 
to-volume ratio which is minimal for spheres. It can 
therefore be concluded that a non-spherical particle of 
the same size as a spherical one has larger values for 
the basic parameters q/m and /3, or, a non-spherical 
particle with the same q/m and /3 values used in the 
simulation is larger than indicated in table 5. Still, 
q/m, and f3 are quantities that depend on material 
properties of the grains and have not been measured 
directly. Therefore, they have to be treated as free 
parameters. 

Another assumption made in the model is the con- 
stancy of q. Due to short-term variations in the solar 
wind, the surface potential and thus q is expected to 
fluctuate. These charge fluctuations have the same 
effect on the grain's motion as short-term variations 



of the ambient magnetic field. As shown above, it is 
a good approximation to use an average strength of 
the Lorentz-force, with the averaging time-scale being 
the Larmor-pcriod of the grain. Furthermore it was 
shown by Zook et al. [1996], who modeled the motion 
of very small Jupiter stream particles through inter- 
planetary space, that the effect of charge fluctuations 
on grain dynamics can be neglected for the grain sizes 
under consideration. 

Finally it can be argued that the solar wind mag- 
netic field is not well represented by a Parker spiral 
during the solar maximum. The physical situation 
during solar maximum is a highly disordered field 
with the polarity distribution showing no clear sep- 
aration of polarities between the northern and south- 
ern heliographic hemisphere [Hoeksema, 1985]. As 
shown above, the motion of grains in the size range 
of 0.1 /im depends only on the average polarity p. In 
the highly disordered configuration of the solar max- 
imum, regions of positive and negative polarity act 
on the grains in a stochastic way, independent of the 
heliographic latitude of the grain's position. There- 
fore p ~ and thus vanishing average Lorentz-force 
for all heliographic latitudes is a good approximation 
of the physical situation during the solar maximum, 
independent of the actual form of the field lines. 

3. Spatial Distributions of Interstellar 
Dust in the Solar System 

The spatial distribution of interstellar grains inside 
the heliosphere shows how the concentration of grains 
of different sizes is reduced or enhanced by solar grav- 
ity, radiation pressure, and the solar wind magnetic 
field. To display the distributions that result from 
the model described in section 2, a coordinate-system 
is chosen such that the x-axis is anti-parallel to the 
initial velocity vector. The z-axis is then chosen to 
be perpendicular to the a;-axis and to lie in the plane 
of the initial velocity vector Voc and the solar rota- 
tion axis vector luq (see figure 1). The y-axis is set 
to be perpendicular to the x- and z-axes such that a 
right-handed coordinate-system is formed. 

The angle k between the downstream-direction and 
the solar rotation axis (as shown in figure 1), has been 
determined by fitting the directional information of 
the Ulysses in-situ measurements. The best fit was 
achieved at k = 98° [Frisch et al, 1999]. 

The result of the simulation is shown for four dif- 
ferent grain radii: 0.7 /im, 0.5 /im, 0.3 /im, 0.2 /im, 
and 0.1 /im. For the dynamical parameters (3 and 



q/m of these grains see table 2. 

The spatial distribution of the 0.7 ^m and 0.5 /um 
grains is shown in figure 2. Because of the small 
charge-to-mass ratio of these grains, the coupling to 
the time-variable solar wind magnetic field is weak, 
and therefore their spatial distribution is nearly con- 
stant and rotationally symmetric about the x-axis 
(downstream direction). The density enhancement 
downstream of the Sun is produced by gravitational 
bending of the grain's trajectories towards the Sun. 
Due to their smaller /3-value, the downstream con- 
centration increase of the 0.5 /xm-grains is less pro- 
nounced. The comparison with the result of the ana- 
lytical calculation [Landgraf and Milller, 1998], which 
assumes q/m = 0, shows that electromagnetic inter- 
actions can indeed be neglected for these particles. 
The spatial distribution of grains with radii of 0.3 /im 
is already slightly affected by their interaction with 
the solar wind magnetic field. 

For even smaller grains the interaction with the so- 
lar wind magnetic field becomes important, as shown 
in figure 3, in which the time-sequence of the den- 
sity distribution of interstellar 0.2 /im-grains in the 
a;-j/-plane (approximately in the ecliptic plane) is dis- 
played. As shown in table 1, the solar maximum 1991 
was the end of a focusing cycle. The density enhance- 
ment in the x-y-plane is the result of the focusing of 
grains to this plane. At the maximum of solar activity 
the average magnetic field is weak, because the cur- 
rent sheet extends to high hcliographic latitudes. In 
this configuration sectors of positive and negative po- 
larity cancel each other efficiently as they pass by the 
dust grains, creating a weak average Lorentz-force. 
Therefore, the region of enhanced grain density moves 
downstream without much perturbation until in 1996 
the defocusing effect becomes important. The void 
region downstream of the Sun is expanded by the de- 
focusing, and the overall density in the x-j/-plane is 
reduced until the end of the defocusing cycle in 2002. 
After that the new focusing cycle starts to increase 
the density in the x-y-plane until the distribution in 
2012 is similar to the one in 1991 when the defocusing 
cycle started. As a consequence of the large /3-value 
of /3 = 1.4, the concentration of 0.2 y^m-grains down- 
stream of the Sun is very low. In the x-^-cut (plane 
approximately perpendicular to the ecliptic that also 
contains the stream vector) of the density distribu- 
tion shown in figure 4 one can see that the density 
distribution is no longer symmetric about the x-axis 
for 0.2 /^m-grains. The first panel in figure 4 also 
shows that because the upstream direction lies close 



to but not in the plane of the solar equator, the dis- 
tribution is asymmetric about the x-y-plane. At the 
end of the focusing cycle in 1991 the grains are con- 
centrated along the plane of the solar equator and the 
density at higher latitudes is reduced. As the defo- 
cusing cycle proceeds, the density at higher latitudes 
increases as more grains are deflected into this region. 
In 2002, at the end of the defocusing cycle the re- 
gions of enhanced density move downstream and the 
new focusing cycle starts again to concentrate grains 
around the plane of the solar equator until in 2012 
the distribution of 1991 is restored. 

The spatial distribution of grains with radii of 
0.1 /xm is dominated by the interaction of these grains 
with the solar wind magnetic field due to their high 
charge-to-mass ratio of q/m =■ 5.31 C kg^ . At the 
end of the focusing solar cycle, the grains are concen- 
trated in a sheet about the plane of the solar equator 
as can be seen in the first panel in figure 6. Because 
the solar wind magnetic field moves radially outward 
with the solar wind, the Lorentz-force on the grains 
is not perfectly perpendicular to the solar equator, 
but also contains a small outward directed compo- 
nent. This force decelerates small grains and causes 
a region of enhanced density upstream of the Sun, 
which is visible as the arc-shaped feature in the first 
panel of figure 5. As the average field strength is low 
during the solar maximum, this feature can propa- 
gate towards the Sun and the concentration around 
the plane of the solar equator is split into a northern 
and southern part by the new defocusing field (1993 
to 2001 in figure 6). As can be seen in the panel for 
1996 in figure 5, the region of high grain density is 
diverted before it reaches the solar vicinity. There- 
fore, the flux of 0.1 /im grains measured in the solar 
vicinity (roughly inside Jupiter's orbit) is expected 
to be reduced during the full solar cycle due to the 
interaction with the solar wind magnetic field. This 
deficiency of small grains has been observed in the 
in-situ data \Grun et at, 1994; Landgraf et at, 1999]. 
After a new focusing cycle starts in 2002, grains are 
again concentrated around the solar equator as can 
be seen in the panels from 2006 to 2012 in figure 6 
until the configuration of 1991 is restored. 

The simulation shows that the flux of small in- 
terstellar grains in the solar vicinity depends on the 
phase of the solar cycle. Since in the interstellar 
medium, small grains (a < 0.2 /im) dominate the 
dust population by number as indicated by extinc- 
tion measurements [Mathis et al., 1977; Kim et al., 
1994], the total flux of interstellar grains in the So- 



lar System should be dominated by smallest grains 
that are not filtered by the solar wind magnetic field. 
Therefore, the total flux of interstellar dust depends 
on the phase of the solar cycle. In the next section 
the flux of interstellar grains measured in-situ with 
the Galileo and Ulysses dust detectors is compared 
with the temporal variation predicted by the model. 

4. Comparison with In-Situ 
Measurements 

The magnitude /isd(ii,^2) of flux of interstellar 
dust measured by the dust detectors on-board Ulysses 
and Galileo in the time-interval [^1,^2] is determined 
by /isd = Nisd{ti,t2)/I{ti,t2), where Nisdih^h) is 
the number of interstellar grains detected in the 
given time- interval and /(ti,t2) is the average sen- 
sitive detector-area, accumulated between ii and ^2- 
I{ti,t2) is given by: 
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where 0rot is the spacecraft rotation angle (Ulysses 
and Galileo are spin-stabilized spacecraft, and there- 
fore the sensitive area for a dust stream from one di- 
rection has to be averaged over the rotation angle), 
A is the sensitive detector-area of the dust detector 
which is a function of time and rotation angle, Vqc is 
the dust velocity at large distances from the Sun, and 
iirci is the relative velocity of the spacecraft and the 
dust as a function of time. By introducing the factor 
I'rci/t'oo, the effect of the flux measured in the space- 
craft frame being enhanced if the spacecraft moves 
upwind and reduced during downwind motion is taken 
into account. Therefore, /isd is the flux in the inertial 
frame at the location of the spacecraft at time t. 

In the following fisd{t) predicted by the model is 
compared with the flux measured with the Ulysses 
and Galileo spacecraft. For the comparison it is as- 
sumed, that the population of interstellar dust grains 
is dominated by grains out of a narrow size-interval 
which can be represented by one of the simulated 
grain-sizes (and thus dynamical parameters /3 and 
q/m). In flgure 7 the model prediction for the flux 
of interstellar grains at the location of the Ulysses 
and Galileo spacecraft is shown as a function of time. 

For the Galileo measurements the simulation pre- 
dicts a constant flux of interstellar grains on the de- 
tector for all but the smallest grain sizes. This is 
because of the short period of time covered by the 
Galileo data. Interstellar grains could be identified 



clearly in the Galileo data only after it left the in- 
ner Solar System in mid-1993 and before it reached 
the Jovian system at the end of 1995. This time in- 
terval is short compared to the duration of the solar 
cycle, therefore the change of the solar cycle did not 
affect the fiux of 0.2 /xm- and 0.3 /im-grains. The 
predicted increase in the fiux of 0.1 fim-giams is due 
to to Galileo's upstream directed motion towards the 
region of enhanced spatial density as can be seen in 
the first three panels of figure 5. 

The prediction for the long-term measurements of 
Ulysses from February 1992 (Jupiter fiyby) to end of 
1997 (latest data) shows that for 0.1 /im-, 0.2 /im, 
and 0.3 ^m-grains the effects of the solar cycle should 
modulate the measured flux. In the Ulysses data, the 
fiux of 0.1 y^m-grains is predicted to be strongly re- 
duced. After 1996 the simulation predicts a decrease 
in the flux of the 0.2 /um- and 0.3 /^m-grains. This is 
due to the defocusing effect of the solar wind magnetic 
field on the trajectories of these grains. The flux of 
larger grains is nearly constant and close to its value 
at large distances, because they couple weakly to the 
solar wind magnetic field. 

Figures 8 and 9 show the comparison of the tempo- 
ral variation of the flux predicted by the simulation 
with the Galileo and Ulysses in-situ measurements, 
respectively. Since the spatial density of grains rep- 
resented by one of the simulated grain sizes at large 
heliocentric distances is not known, the normalization 
of the spatial distribution is a free parameter that can 
be scaled by an overall factor such that the square of 
the difference of prediction and measurement is min- 
imized. If grains with dynamical parameters similar 
to the simulated 0.1 y^m-, 0.2 fim-, or 0.3 /im-grains 
dominate the interstellar dust population measured 
by Galileo and Ulysses, the predicted temporal varia- 
tion of the flux of the corresponding grain size should 
match the measurements. 

Galileo has measured an increasing flux of inter- 
stellar grains on its way to Jupiter. According to the 
simulation this is best explained with grains that in- 
teract strongly with the solar wind magnetic held and 
are thus deficient in the inner Solar System. Because 
Galileo was moving upwind to larger heliocentric dis- 
tances, it detected an increasing fiux. The same de- 
ficiency of grains in the inner Solar System can be 
created by the effect radiation pressure, but to ex- 
plain an increase of fiux outside 4 AU observed by 
Galileo, the grains must have very high (3 values of 
2.5 and larger, which is not supported by the Mie- 
calculations. Compared to the duration of the solar 



cycle, Galileo spent a short time in the interplanetary 
space of the outer Solar System. Therefore, the com- 
parison of the Galileo measurements with the results 
of the simulation is rather ambiguous and does not 
give any information about the effect of the chang- 
ing solar cycle on the distribution of interstellar dust 
grains in the Solar System. 

Since the Ulysses measurements cover nearly 6 
years, the effect of the changing solar cycle can be ob- 
served in the Ulysses data. The measured flux stayed 
constant after Ulysses left the ecliptic plane in Febru- 
ary 1992 until mid-1996. After that the flux of inter- 
stellar grains decreased until it reached its minimal 
value of (4± 1) • 10"'"^ m^^ g-i at the end of 1997. For 
0.1 lira grains the simulation does not predict a con- 
stant flux at the Ulysses location between 1992 and 
1996. Furthermore, the flux of the simulated 0.1 /itm- 
grains decreases by an order of magnitude within one 
year after mid-1996. This is also not observed. There- 
fore, the in-situ measurements as well as the simula- 
tion indicates that grains with charge-to-mass ratios 
of 5 C kg~^ or larger do not contribute significantly 
to the flux of interstellar grains in the Solar System 
within Jupiter's orbit. For 0.2 /xm- and 0.3 /im- 
grains (charge-to-mass ratios between 0.5 C kg^ 
and 1.4 C kg~ ) the simulation predicts a flux varia- 
tion that is in good agreement with the observation, 
though long-term measurements are necessary to con- 
flrm the agreement. Interstellar grains with lower 
charge-to-mass ratios do not show a decrease in flux 
(see right panel of flgure 7) . This explains why bigger 
grains dominate the mass distribution of interstellar 
grains measured by Ulysses after mid-1996 compared 
to the mass distribution measured between 1992 and 
mid-1996 [Landgraf et al, 1999]. 

5. Conclusion 

The model described in section 2 explains the ob- 
served time variability of the flux and grain mass dis- 
tribution of interstellar dust in the Solar System as 
due to electromagnetic interaction of the grains with 
the changing solar wind magnetic field. It also shows 
that very small grains with large charge-to-mass ra- 
tios are repelled from the Sun by the average elec- 
tromagnetic effects. Since the grain size distribution 
in the LIC is believed to be very steep, the most 
abundant interstellar grains have charge-to-mass ra- 
tios just small enough (or, equivalently, sizes large 
enough) for their inertia to overcome the electromag- 
netic repulsion. In this marginal case the electromag- 



netic perturbation is still strong and defines the spa- 
tial distribution of the grains. The knowledge of the 
spatial distribution of interstellar dust grains can be 
used to predict in-situ measurements by spacecraft 
(e.g. Cassini and Stardust) [Landgraf and Muller, 
1998] depending on their location in the Solar Sys- 
tem. 

The simulation has furthermore shown that the 
spatial distribution of small interstellar dust grains 
deviates strongly from the distribution of bound, in- 
terplanetary grains. This asymmetry as well as the 
time variability of the spatial distribution can be used 
to observe the infrared emission of interstellar dust in 
the Solar System by telescopic means and distinguish 
it from the radiation emitted by interplanetary grains. 
Because of the low average spatial density (about 
10""'^^ cm~^ to 10~^^ cm^"^) of interstellar grains, an 
instrument of very high sensitivity is needed for such 
an observation. 

Further long-term in-situ measurements are needed 
to confirm the time variation of the interstellar dust 
flux discovered in the Ulysses data. The simulation 
has shown that even small interstellar dust grains 
reach heliocentric distances of 1 AU. Therefore, an 
Earth orbiting satellite can be used to measure or 
even collect dust grains from outside the Solar Sys- 
tem. 
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Figure 1. Coordinate-system in which the spatial distribution of interstellar dust in the heliosphere are displayed. The 
a;-j/-plane is close (< 10° off) to the cliptic. A tilt angle k of 90° between the downstream-direction and the solar rotation 
axis is compatible with the determination of the stream-direction, but the best fit to the in-situ data is achieved at k = 98° 
[Frisch et al, 1999]. 
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Figure 2. The spatial distribution of 0.7 fiui- (upper two panels) and 0.5 fim- (lower two panels) grains in a 20 AU x 20 AU 
area in the x-y-plane. The linear color scale represents the spatial grain density in units of the density at infinity. The two 
panels on the left have been calculated using the numerical simulation of particle trajectories and the two panels on the right 
show the result of a analytical calculation for q/m = but the same /3- values (see table 2). The numerical simulation is in 
good agreement with the results from the analytical calculation, the deviations result from the limited number of particles 
and the limited spatial resolution of the simulation. As expected, gravitational focusing is weaker for 0.5 /im-grains due to 
their higher /3- value. 
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Figure 3. Distribution of simulated interstellar 0.2 /xm grains in the a::-2y-plane during the solar cycle between 1991 and 
2012. The time sequence starts at the upper left panel and ends at the lower right. Time steps are about 2.5 years. Plots 
are 80 AU x 80 AU 
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Figure 4. Distribution of simulated interstellar 0.2 /xm grains in the x-z-plsme during the solar cycle between 1991 and 
2012. The time sequence starts at the upper left panel and ends at the lower right. Time steps are about 2.5 years. Plots 
are 80 AU x 80 AU. 
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Figure 5. Distribution of simulated interstellar 0.1 /im grains in the a::-2y-plane during the solar cycle between 1991 and 
2012. The time sequence starts at the upper left panel and ends at the lower right. Time steps are about 2.5 years. Plots 
are 80 AU x 80 AU. 
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Figure 6. Distribution of simulated interstellar 0.1 nm grains in the x-z-plsme during the solar cycle between 1991 and 
2012. The time sequence starts at the upper left panel and ends at the lower right. Time steps are about 2.5 years. Plots 
are 80 AU x 80 AU. 
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Figure 7. The flux of interstellar dust in the inertial frame at the location of Galileo (left panel) and Ulysses (right panel) 
as a function of time as predicted by the simulation. The magnitude of the flux is normalized to its value at large distances 
from the Sun. The radius a of the simulated grains increases from the top to the bottom. Major mission events are indicated 
by dotted vertical lines. In both cases, the flux of small grains varies stronger than the flux of big grains due to their coupling 
to the solar wind magnetic field. The flux of 0.1 /im-grains strongly reduced at Ulysses locations. 
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Figure 8. Comparison of the temporal variation of the flux in the inertial frame at Gahleo's location predicted by the 
simulation (solid line for 0.1 /im-grans, dashed line for 0.2 /im-grans, and dash-dotted linefor 0.3 /im-grans) with Galileo 
measurements (crosses). Due to the short time interval covered by the measurements, no significant changes of the flux due 
to the solar cycle are expected. 
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Figure 9. Comparison of the temporal variation of the flux in the inertial frame at Ulysses' location predicted by the 
simulation (solid line for 0.1 /im-grans, dashed line for 0.2 /im-grans, and dash-dotted linefor 0.3 /im-grans) with Ulysses 
measurements (crosses) . Because impacts of interstellar grains could not be distinguished clearly from impacts of interplan- 
etary grains around perihelion passage, there is a gap in the measured flux between January and May 1995. The measured 
decrease in flux after 1996 can be due to the onset of the defocusing cycle in 1991 as indicated by the results of the simulation. 
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Table 1. Average polarity phases of the solar wind mag- 
netic field in the 22-year solar cycle parameterized by the 
current sheet tilt angle e. If the north/south polarity con- 
figuration is -I-/-, the magnetic north pole is the helio- 
graphic north pole, and for -/-f the magnetic north pole 
is the heliographic south pole. 
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Table 2. Dynamic parameters of the simulated grain sizes. The 
mass is determined by assuming the grains to be homogeneous 
spheres with a density of 2.5 g cm^"^. The radiation pressure parame- 
ter /3 is taken from Mie-calculations for spheres made of astronomical 
silicates [Gustafson, 1994; Draine and Lee, 1984], and the charge is 
calculated for a sphere with a surface potential oi U — +5 V. 
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